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The statistical mechanics of single polymer knots is studied using Monte 
Carlo simulations. The polymers are considered on a cubic lattice and their 
conformations are randomly changed with the help of pivot transforma- 
tions. After each transformation, it is checked if the topology of the knot is 
preserved by means of a method called pivot algorithm and excluded area 
(in short PAEA) and described in a previous publication of the authors. 
As an application of this method the specific energy, the radius of gyration 
and heat capacity of a few types of knots are computed. The case of attrac- 
tive short-range forces is investigated. The sampling of the energy states is 
performed by means of the Wang-Landau algorithm. The obtained results 
show that the specific energy and heat capacity increase with increasing 
knot complexity. 

PACS numbers: 61.41. +e, 65.60.+a, 61.43.Bn 

1. Introduction 

Polymer knots are researched in connection with several applications, 
mainly in biology and biochemistry [BEIEIHIEIISIEIEIEIIES. In this 
work we investigate the statistical mechanics of a single polymer knot by 
computing its specific energy, heat capacity and gyration radius at differ- 
ent temperatures. One of the main problems in studies of this kind is how 
to preserve the polymer topology during the thermal fluctuations. To this 
purpose, several methods have been developed [HI [12l [131 El E3 13 HZ|- 
Most of them are based on the Alexander or Jones polynomials, which are 
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rather powerful topological invariants able to distinguish with a very high 
degree of accuracy the different topological configurations. The main draw- 
back of these polynomials is that their calculation is time consuming from 
the computational point of view. In a recent work [18] , a strategy based on 
an excluded area method called PAEA has been proposed in order to cir- 
cumvent these difficulties. A similar idea has been presented in |19j . but in 
that case, instead of using arbitrary pivot transformations to produce ran- 
dom knot configurations, a set of topology preserving pull moves is adopted. 
Within our method, instead, any transformation is allowed. In this way the 
equilibration of the system is faster and the access to all possible conforma- 
tions is easier. Those transformations that lead to a change of topology are 
automatically detected by the PAEA algorithm and discarded. 

According to |18j . one starts from a seed knot, for instance those given 
in ref. [20]. Next, the knot is changed at a randomly chosen set of segments 
by applying the pivot algorithm [21J. After each pivot move, it is easy to 
realize that the difference between the old and new configurations, obtained 
by canceling the segments that have been unaffected by the transformation, 
consists of a closed loop. Around this loop, an arbitrary surface is stretched, 
whose boundary is the closed loop itself. The criterion to reject changes that 
destroy the topology of the knot is the presence or not of lines of the old 
knot that cross such surface. If these lines are present, the trial pivot move 
is rejected, otherwise is accepted. This combination of pivot algorithm and 
excluded area (PAEA) provides an efficient and fast way to preserve the 
topology that can be applied to any knot configuration, independently of 
its complexity. For pivot moves involving a small number of segments the 
method becomes exact. This technique may be employed in the study of 
the thermal and mechanical properties of polymer knots as it has been done 
in Ref. [18] . In this article we will extend that work by studying the case 
of an attractive short-range potential, which is nothing else but a rough ap- 
proximation of the Lennard-Jones potential on the lattice. Moreover, with 
respect to [18] . we do not limit ourselves to the computation of the internal 
energy and heat capacity, but we consider also the gyration radius of the 
knot. The cases of the unknot, trefoil and 5i topologies is investigated. The 
sampling of the canonical ensemble is achieved by using the Wang-Landau 
algorithm |22j at different temperatures. 



2. Sampling and calculation method 

The way of generating random knot transformations with the help of 
pivot transformations and the PAEA method needed to prevent topology 
changes after these transformations have been already extensively described 
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in [21] and p2] respectively. We refer the interested reader to those publica- 
tions for more details. In this Section we concentrate on the Wang-Landau 
(WL) method applied to polymer knots. See also [23] for applications of 
this method to linear polymer chains and rings. 

The WL algorithm can be regarded as a self-adjusting procedure for 
obtaining the density of states 

J dX6(Ei - E(X)) = (1) 

where X is a microstate of the system under consideration. We suppose 
here that the energy values are discrete, so that they can be labeled by 
indexes If the $j's are known, then the partition function can 

be constructed: Z = Y,i e ~ PEi ®^\ 

As well, it is possible to derive in an 
easy way the averages of any quantity that can be expressed in terms of 
the momenta of the energy (E 1 )(J3) = Y J i E \ e ~ fiE ^il z ^ Z = 1,2,.... For 
example, the heat capacity is given by: 

C(T) = p\(E 2 )(P)-((Emf) (2) 

According to the Wang-Landau method, the density of states is constructed 
with successive approximations. First, the would be density of states g(Ej) 
is set to be equal for all E^s by putting g(Ei) = 1. After that, a Markov 
chain of microstates X%,X2, ... is generated. The probability of transition 
from a state Xi of energy Ei to a state with energy is 

p(i^i') = mm[l,^\] (3) 
9{Ev) 

If p{i —> i') > 1, the state i! is automatically accepted. If instead p(i — > i') < 
1, then a random number < rj < 1 is generated and the state is accepted 
only if p(i —> i') > r]. Once an energy state Ei is visited, its corresponding 
would be density of states is updated by multiplying it by a modification 
factor /, i. e. 

g(Ei) -> fg(Ei) (4) 

where / > 1. Moreover, the energy histogram H(Ei) is updated by perform- 
ing the replacement H(Ei) —> H(Ej) + 1 |22j. In this way, an energy state 
occurring N{Ei) times during the sampling will have a would be density of 
states g(Ei) = f N ( E d and the transition probability Q to that state will be 
suppressed by the factor f~ N ( Ei \ When the energy histogram H(E{) be- 
comes flat, then g(E{) converges to the density of states To show that, 

We have put here f3 — T , where /3 is the usual Boltzmann factor in thermodynamic 
units in which the Boltzmann constant is equal to 1. 
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let's consider the probability V{Ej) of obtaining a microstate X with energy 
E(X) = E{. This probability must be equal to the probability of generating 
the state X times the probability of acceptance of E{ introduced by the 
Wang-Landau algorithm, which is proportional to f~ N ( E i), I n formulas: 

P(Ei) oc r Nm ^ (5) 

The last factor in the above equation is due to the fact that the probability 

of obtaining a microstate X with energy E{ is given by jJ d j X d xf^ E ^^(x ) ) = 

•^2 ^. and the denominator ^ <3?j is an irrelevant constant. When the energy 

histogram H{Ei) becomes flat, this means that the probability V{Ei) is the 
same for all states E%. In other words: V{Ei) = a for every i, where a is a 
constant. Thus, from Eq. (|5| we obtain: 

f N W = giEt) oc <5>i (6) 

Actually, if / is too big, the statistical errors on the g(Ei) y 8 may grow 
large and the above equation is satisfied very roughly. On the other side, 
if / is too small, it is necessary an enormous number of microstates during 
the sampling in order to derive the g(Ei)'s. For this reason, in the Wang- 
Landau method the density of states is computed perturbatively. In the 
next step, one takes the g(E{) evaluated with the modification factor / 
as a starting point and generates another Markov chain of microstates. 
The energy histogram H{Ei) is reset to and all the previously explained 
procedure is repeated for the new microstates apart from the fact that in 
Eq. Q / is replaced by y/J. When the energy histogram becomes flat, the 
second approximation of the g{Ei) J s is obtained. In the next approximations 
successive square roots of / are entering in the algorithm until we arrive at 
a step n such that yff = ffi na \ ~ exp(10~ 8 ) [22J. The initial parameter / 
is chosen in such a way that the simulations will not take too much time 
and the statistical errors on the g(Ei) will not be too large. 

In the present article the states are distinguished by the number m 
of closest contacts between the monomers, where m takes positive inte- 
ger values. The meaning of contact in the present context is explained in 
Refs. [181 E31 - Short-range attractive forces are studied, so that the en- 
ergy values are given by E m = me, where e is the contact energy of two 
unbonded monomers, which is negative in the attractive case. Since the 
number of samples is huge for polymer systems, it is more convenient to 
consider the logarithm of the density of states O m = In g(E m ). In this way 
the transition probability is expressed as: 

p(m — > m) = min[l, exp(ri m — f2 m ')] (7) 



appb2 printed on January 1, 2013 



5 




Fig. 1. Specific energy, specific heat capacity and gyration radius for the unknot 
and the trefoil as functions of the temperature. Polymers with lengths L = 70, 90 
and 120 are considered, (a) Plot of the specific energy (in e-units); (b) Plot of the 
specific heat capacity (in e-units); (c) Plot of the radius of gyration. 



and the modification factor becomes ln(/). After a state E m is visited, the 
corresponding energy histogram should be updated by H(E m ) — > H(E m ) + l 
and the density of state is modified by Q m — > Q, m + ln(f). 



3. Thermal properties of polymer knots 

In this Section the thermal properties of a few types of polymer knots 
are studied. In particular, the specific internal energy of the polymer per 
unit of length (E)(f3)/L, the heat capacity and the radius of gyration are 
computed in the case of the unknot 0i, the trefoil 3i and the knot 5j.. 

The gyration radius is not directly related to the moments of the energy 
as mentioned in the previous Section. However, this quantity may be com- 
puted by noticing |24j that the mean square radius of gyration (Rq)((3) can 
be written as follows: 

{Rg){P) ~ V r-Pmer, 

Here (-R^m = (j^ J2i j = i((Ri—Rj) 2 )m denotes the average of the gyration 
radius computed over states with m contacts. Moreover, Rj is the position 
vector of I-th segment and L is the length of the polymer. 

In Fig. [T] are displayed the results for the unknot and the trefoil. It 
is found that the growth of the specific energies (E(j3))/L in Fig. [ija) is 
characterised by three regions. At very low temperatures the energy growth 
is practically zero because the temperatures are too low to allow contacts 
between the monomers. When the energy is enough to excite more states, 
the specific energy grows rapidly as a function of the temperature until 
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Fig. 2. Specific energy and heat capacities as functions of the temperature for the 
unknot, the trefoil and the knot 5i with length L = 90. (a) Plot of the specific 
energy (in e-units); (b) Plot of the specific heat capacity (in e-units); (c) Plot of 
the gyration radius. 

saturation is reached and the energy increase becomes moderate. The fast 
increasing energy region causes a peak of the heat capacity as shown in 
Fig. [ljb) as it has been explained in details in [18]. Concerning the topo- 
logical effects, we observe that both the energy and the heat capacity grow 
with growing knot complexity, as it turns out from Fig. [2] by comparing 
the plots for various knots with the same length L = 90. Fig. [TJc) shows 
that in the attractive energy case the mean square gyration radius (Rq)((3) 
grows with increasing temperatures. This is an expected behavior. As a 
matter of fact, in the case of attractive forces the polymer turns out to be 
in a crumpled conformation at very low temperatures with many contacts 
in order to minimize the energy. On the contrary, at high temperatures the 
energy of the thermal fluctuations is large with respect to e, so that the 
attractive forces become negligible. Thus, the number of contacts m will be 
in the average smaller at higher temperatures than at lower energies. As it 
is intuitive, a smaller number of contacts m corresponds to a larger volume 
occupied by the knot, which causes the observed increase of the radius of 
gyration with growing temperatures. 

4. Conclusions 

We have studied the thermal properties of a few types of polymer knots 
using the PAEA algorithm of Ref. [18]. The number of polymer segments 
affected by the pivot moves has been limited to four. In this case the 
PAEA method is able to preserve the topology of the knot exactly |18j . 
The results of |18j . which took into account only short-range repulsive forces, 
have been extended to attractive forces and the calculation of the gyration 
radius has been added. The thermal properties of the unknot, the knot 3i 
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and 5i have been analysed with the help of Monte Carlo simulations based 
on the Wang-Landau algorithm and the pivot method. A brief account 
about the Wang-Landau method has been provided. The results, including 
those with the new topology 5i, confirm those of Ref. [18] even when the 
interactions are attractive. In particular, the presence of the three regimes 
of growth of the specific energy mentioned in the previous Section has been 
observed. Moreover, the role of the topological effects, which make both the 
energy and the heat capacity increase with increasing knot complexity, is 
confirmed. The data coming from the calculation of the gyration radius give 
a measure of the size of the polymer knot at different temperatures. When 
the temperature is low, the number of contacts is at its maximum and the 
gyration radius is at its minimum. The size of the polymer points out to a 
possible crumpled conformation. At high temperatures, the influence of the 
attractive forces becomes negligible and the gyration radius attains slowly its 
maximum. Even if we limited ourselves to simple short-range interactions, 
there are no obstacles to extend our procedure to more realistic polymer 
systems. 
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